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5.4 POLYNOMIAL WIGNER-VILLE DISTRIBUTIONS 0 

5.4.1 Polynomial FM Signals 

A complex signal can be written in the form 

z{t) = a(t)e?W\ (5.4.1) 

where the amplitude aft) and instantaneous phase <j)ft) are real. Let us define 

Ut) = ® • (5.4.2) 

This f z (t) is the instantaneous frequency (IF) of z(t) if z(t) is analytic. In this 
Article we simply use Eqs. (5.4.1) and (5.4.2) as definitions of zft) and f z (t), without 
assuming that z(t) is analytic. Results concerning f z (t) will also be valid for the IF 
[usually written fi(t)] when z(t) is analytic. 

The factor aft) allows amplitude modulation, phase inversion and time limiting. 
If cj>ft) is a polynomial function of degree p, so that f z (t ) is a polynomial of degree 
p — 1, then z(t) is a polynomial-phase or polynomial FM signal. If p > 2, then 
f z ft) is nonlinear, so that z(t) is an example of a nonlinear FM signal. 

Such nonlinear FM signals occur both in nature and in man-made applica- 
tions [1]. For example, the sonar systems of some bats use hyperbolic and quadratic 
FM signals for echo-location. Some radar systems use quadratic FM pulse compres- 
sion signals. Earthquakes and underground nuclear tests may generate nonlinear 
FM seismic signals in some long-propagation modes. The altitude and speed of an 
aircraft may be estimated from the nonlinear IF of the engine noise reaching the 
ground. Nonlinear FM signals also appear in communications, astronomy, teleme- 
try and other disciplines. As these examples suggest, the problem of estimating the 
IF of a nonlinear FM signal is of some practical importance. 

For a deterministic linear FM signal, the Wigner-Ville distribution (WVD) gives 
an unbiased estimate of the IF. To obtain the same property with higher-order 
polynomial FM signals, an extension of the WVD called the polynomial Wigner- 
Ville distribution (PWVD) was defined [2,3]. If the instantaneous phase (f>ft) is 
a polynomial of degree not exceeding p, then the IF estimate based on a PWVD of 
order p is unbiased [4] . 

5.4.2 Principles of Formulation of Polynomial WVDs 

We seek a Time-Frequency Distribution of the form 

P z (t,f) = JF f {R z (t,r)} (5.4.3) 
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where R z (t,r ), called the signal kernel , 1 somehow depends on z(t). If a(t) = 1, 
then, for the best possible representation of the IF law, we would like P z (t,f) to 
be equal to S(f — f z (t)). Making this substitution in Eq. (5.4.3) and taking inverse 
FTs, we find 

R z (t, t) = e j2nf ^ T = e>*'W T . (5.4.4) 

So, while z(t) may have a variable frequency with respect to t, we want R z (t,r) to 
have a constant frequency w.r.t. r, namely f z {t). That is, we want the signal kernel 
to dechirp the signal, yielding a sinusoid of constant frequency, for which the FT 
is the optimal representation and gives a delta function at f z (t). 

To estimate 4>'(t) in Eq. (5.4.4), we use a central Unite-difference (CFD) approx- 
imation of the form 

1 9 

~ <P(t) = - bi[c/)(t+ciT) - (p{t qt )] (5.4.5) 

7 i=i 

and choose the dimensionless 2 real coefficients 6; and Ci so that the approximation 
is exact if <j>(t) is a polynomial of degree p. Let that polynomial be 

p 

(j)(t) = yv (5.4.6) 

i= 0 

so that 

p 

(j) (t) = . (5.4.7) 

»= l 

A polynomial of degree p remains a polynomial of degree p if the argument is 
shifted and scaled. If we shift the time origin so that t = 0 in Eqs. (5.4.5) to (5.4.7), 
we see that if <p(t) contains only even-power terms, both and (j>'{t ) are zero, 
so that the even-power terms do not introduce any errors into the estimate. So we 
may assume that p is even and consider only the odd-power terms in <j)(t). There 
are p/2 such terms, hence p/2 degrees of freedom in the coefficients a, , suggesting 
that the CFD estimate can be made exact by using only p/2 finite differences in 
Eq. (5.4.5), with uniform sampling intervals. But we shall retain the generality of 
the sampling instants because we can do so without further algebraic complexity, 
and because the extra degrees of freedom turn out to be useful. With q = p/2, 
Eq. (5.4.5) becomes 



0'(t) = 



1 

r 



p/2 

Y, bi[<j)(t + qt ) - 4>{t - C t T )] 
1=1 



(5.4.8) 



1 This notation is consistent with the convention that R z (t,r) [Article 2.7] is a generalization 
of K z (t , r) [Section 2.1.2]. But here, as we shall see, the generalization is in a new direction. 

2 We could allow ci to have the dimensions of time and dispense with the symbol r; however, 
retaining r will emphasize the correspondence between the PWVD and the ordinary WVD. 
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where the use of ft instead of ft asserts the exactness of the estimate. Substituting 
this into Eq. (5.4.4), and renaming the signal kernel as R^\t, r) to acknowledge the 
dependence on p, we obtain 

p / 2 b[ 

R^(t, t) = H [ e ^(‘+ clT ) e~ m ~ clT) ' . (5.4.9) 

i=i 

Then, substituting from Eq. (5.4.1) with a(t) = 1, we find 

p/2 

R^\t, r) = U [. z(t+ciT ) z* (t — ciT )] hi . (5.4.10) 

i=i 

The resulting TFD, denoted by W z p \t, /) and given by Eq. (5.4.3) as 

WM(t,f) = JF f {R<?\t,T)} , (5.4.11) 

is called a polynomial Wigner distribution (or polynomial WD) of order p. 

Thus we arrive at a general definition: A polynomial WD of order p of the 
signal z(t) is a function W z p \t, f) given by Eqs. (5-4.10) and (5-4-11), such that 
the coefficients bi and ci satisfy Eq. (5-4-8) when is a polynomial of degree 
not exceeding p. In the special case in which z(t) is analytic, the polynomial WD 
becomes the polynomial Wigner-Ville distribution (PWVD). 

If we put p= 2, ci = 5 and bi = 1 , then R^ (f, r) reduces to z(t +^ ) z*(t — |) , 
which is the instantaneous autocorrelation function (IAF), denoted by K z (t, r). We 
might therefore describe Ri p \t,r) as a polynomial IAF or higher-order IAF. 

Eq. (5.4.10) shows that for any z(t), the polynomial IAF is Hermitian in r. It 
follows that the polynomial WD is real. 

5.4.3 IF Estimates with Zero Deterministic Bias 

When the IF is a polynomial of degree not exceeding p— 1, the PWVD of order p 
gives an unbiased estimate of the IF law, as is shown by the following result. 

Theorem 5.4.1: If z(t) and f z (t) are given by Eqs. (5-4-1) and (5-4-2), where 
<f>[t) is a polynomial of degree not exceeding p, and ifW z P \t,f) satisfies the general 
definition of a pth- order polynomial WD of z(t), then wi P \t, f) is symmetrical in 
f about f = f z (t). 

Proof: Substituting Eq. (5.4.1) into Eq. (5.4.10) and simplifying, we find 

p/2 

R { z P \t,r) = R(f\t,r) exp (j^2 bfftt+cir) - ftt-cir)]^ (5.4.12) 

V j = i 

where Rftft, r), as the notation suggests, is the polynomial IAF for a(t), given by 

p/2 

Ri P \t,r) = [a(t+cir) a(t-ciT)] bl . 
l=i 



(5.4.13) 




188 



Chapter 5: Design of Time-Frequency Distributions 



Because f(t) is a polynomial function of degree not exceeding p, and because W z p \t, f) is 
a pth-order polynomial WD, Eq. (5.4.8) is applicable, so that Eq. (5.4.12) becomes 

R ( z p \t,r) = R (p \t,T)e j ' l ’' (t)T = R {p) (f,r) e j2nMt)T . (5.4.14) 

Taking Fourier transforms (r — + /), we find 

Wi p \t,f) = (5.4.15) 

where Wa P \t,f) is the corresponding polynomial WD of a(t): 

W^ p \t,f) = J^{i#>(t,T)}. (5.4.16) 

From Eq. (5.4.13) we see that i?i P ^(t, r) is real and even in r. Hence, from Eq. (5.4.16), 
Wa P \t, f ) is real and even in /. Then, from Eq. (5.4.15), W z p \t, /) is real and symmetrical 
in / about f = fz(t). ■ 

Because a symmetrical distribution is symmetrical about its first moment, The- 

orem 5.4.1 has the following corollary: If z{t) has polynoviial phase of degree not 
exceeding p, and if W~ p \t, f) is a pth- order polynomial WD of z(t), then the first 
moment of W z p \t, f) w.r.t. f is equal to f z {t). 

Being unbiased for higher-degree polynomial FM signals, PWVDs can solve 
problems that quadratic TFDs cannot [2], PWVDs also give optimal frequency 
resolution in the sense that they are FTs of maximal-length polynomial IAFs derived 
from the full-length signal. 



5.4.4 Calculation of Coefficients 

Applying Eq. (5.4.6) in Eq. (5.4.8) gives 



m 



, P/2 



h [(f+c/r) 1 - (t-cir)\ 



1=1 L i= 0 



(5.4.17) 



Eqs. (5.4.7) and (5.4.17) give two expressions for <j>'(t). Equating these expressions, 
and shifting and scaling the time variable so that t = 0 and r = 1, we obtain 



a i 



p/2 

E 

i=i 



p 


I p/ 2 1 


v - 1 


bij^ai W ~ (— ci)*] 


II 

to 

M 


^ ^ Q 


_ i= 0 


i i=i 


2=1,3,... 



h. 



(5.4.18) 



In the left-hand and right-hand expressions of this equation, the coefficient of ai is 
zero for all even values of i, justifying the decision to consider only odd values of i. 
Equating coefficients of ai gives 



1 

2 



p/2 

E c ^ ■ 

i=i 



(5.4.19) 
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and equating coefficients of 03 , a 5 , a p _i gives 



p/2 



o = 5 > f &i; * = 3,5, . . . ,p - 1. 



i=i 



two 


equations 


can 


be written in 


Cl 


C2 


C 3 


Cp /2 


ci 3 


c 2 3 


C3 3 


Cp /2 


Cl 5 


C2 5 


C3 5 


Cp /2 


Ci p - 


- 1 C2?- 1 


C3 P 


_1 ••• Cp/2' 



3 

5 

p-i 





' h - 




■ 1/2 - 




&2 




0 




bz 


= 


0 




- bp/2 . 




0 



(5.4.20) 



(5.4.21) 



and solved algebraically or numerically. 

Eq. (5.4.21) is underdetermined: p/2 equations in p unknowns leave p/2 degrees 
of freedom. Various ways of exploiting the degrees of freedom give rise to various 
forms of the PWVD. Here we mention two forms described in [2, p. 217]. 

If we decide that we want uniform sampling intervals, we choose the values of c; 
and accept the resulting values of bi , which in general turn out to be fractions. The 
resulting form of the PWVD is called the fractional- powers form or Form I, of 
which an example is given in [3, p. 550]. The need to compute fractional powers is 
a cause of inefficiency. 

Alternatively, if we decide that the polynomial IAF must contain only positive 
integer powers, we choose the values of 6/ and accept the resulting values of q , 
which in general give non-uniform sampling intervals. The resulting form of the 
PWVD is called the integer-powers form or Form II. In a discrete-time imple- 
mentation, non-uniform sampling requires interpolation, but this is more efficient 
than computing non-integer powers. The number of interpolations required can be 
reduced by time-scaling [2, p. 218]. 

(There is also a “complex-time form” or “Form III” based on an analytic extension 
of the signal with a complex time argument; this is described in [5], with further 
details in [6].) 

Notice that given one solution to Eq. (5.4.21), we can obtain another solution by 
changing the sign of any c; and of the corresponding 6;. Thus, from one solution, 
we can always find another solution in which all the indices 6/ are positive, ensuring 
that the polynomial IAF does not become unbounded as the signal approaches zero. 

5.4.5 Examples 

In the trivial case for which p = 2 (quadratic phase, linear FM), Eq. (5.4.21) reduces 
to C161 = 1/2. If we fix the sampling points by choosing ci = 1/2, that leaves 
6i=l. Alternatively, if we choose bi = 1 for unit powers, that leaves Ci = 1/2 . 
Thus we have a degenerate case in which Forms I and II overlap. Substituting for 
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p , b\ and c\ in Eq. (5.4.8), we obtain the simplest possible CFD approximation: 

0'(*) = ^ [*(*+§)- *(*-§)]■ (5-4.22) 

The same substitutions in Eq. (5.4.10) yield the polynomial IAF 

R?\t,T) = z{t+l)z*{t-l). (5.4.23) 



This is just the ordinary IAF K z (t,r), which when substituted into Eq. (5.4.11) 
yields the ordinary (quadratic) Wigner distribution. So, for p = 2, Theorem 5.4.1 
says that the Wigner distribution is symmetrical about f z (t) if z(t) is a quadratic- 
phase (i.e. linear FM) signal. This confirms that the WVD gives an unbiased esti- 
mate of the IF for deterministic linear FM signals. 

In the case for which p = 4 (quartic phase, cubic FM), Eq. (5.4.21) reduces to 
the 2x2 system 



Ci&i + C 2 b 2 = 1/2 
Ci 3 &i + c 2 3 b 2 = 0 



(5.4.24) 



Form I: If we take Ci = 1/4 and c 2 = —1/2, Eqs. (5.4.24) become a linear system 
with solutions bi = 8/3 , b 2 = 1/3 . Substituting these values into Eq. (5.4.10) gives 






(5.4.25) 



where the superscript “(4)” indicates order 4. 

Form II: If we take bi = 2 and 6 2 = 1 , Eqs. (5.4.24) become a nonlinear system 
with solutions 



Cl 



2 ( 2 - 2 1/3 ) 



0.6756 ; c 2 = -2 1/3 ci « -0.8512 . 



Substituting for b\ and b 2 in Eq. (5.4.10) gives 

R^(t,r) = [z(t+cir) z*(t-cir)} 2 z(t+c 2 r) z* {t-c 2 r). 



(5.4.26) 



(5.4.27) 



where Ci and c 2 are as in Eqs. (5.4.26). This is one of an infinite number of Form-II 
solutions. More details on the design procedure are given in [1,2,7] and in Article 5.5 
(next). 

Fig. 5.4.1 illustrates the effectiveness of the PWVD in suppressing artifacts 
caused by errors in the CFD estimate of the IF. The IF of the test signal is 



fz(t) = f c + f d C 5 (t/T) ; 0 <t<T (5.4.28) 



where f c is the center frequency, /d is the minimax frequency deviation, and C$ is 
the Chebyshev polynomial of degree 5. Trace (a) shows the ordinary Wigner- Ville 
distribution (WVD), which clearly cannot handle the nonlinear IF law. The Form-II 
PWVD shown in trace (b) was computed by the “ TFSA" software toolbox, developed 
in-house by the Signal Processing Research Centre, QUT. It is of 6th order, so that 
its signal kernel exactly dechirps the 5tli-degree IF law. The superiority of the 
PWVD is evident, as is its symmetry about the IF law. 
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Fs=1Hz N=I29 Fs=1Hz N=129 




Fig. 5.4.1: Time-frequency representations of a 5th-degree Chebyshev polynomial FM signal 

[Eq. (5.4.28)] with center frequency f c =0.25 Hz, minimax frequency deviation /d = 0.1 Hz, dura- 
tion T = 128 seconds: (a) WVD; (b) 6th-order Form-1 1 PWVD. Both TFDs are unwindowed. Each 
plot shows time vertically (range 0 to 128 s; resolution 2 s) and frequency horizontally (range 0 to 
0.5 Hz), and has the time trace at the left and the magnitude spectrum at the bottom. 



5.4.6 Multicomponent Signals and Polynomial TFDs 

The use of TFDs for multicomponent signal analysis requires a reliable method of 
suppressing cross-terms (see Section 3.1.2 and Articles 4.2 and 5.2). 

For quadratic TFDs, the mechanisms of generation and suppression of cross- 
terms are well understood. In the WVD, each cross-term is generated midway 
between the interacting components and alternates at a rate proportional to the 
separation between the components (see Article 4.2). Other quadratic TFDs may 
be obtained from the WVD by 2D low-pass filtering, which attenuates the cross- 
terms because of their alternating (high-pass) character. 

For higher-order polynomial TFDs, there is the added difficulty that cross-terms 
do not necessarily alternate, so that it may not be possible to suppress cross-terms 
entirely by convolving them with a simple smoothing function in the (t, /) plane. 
However, if polynomial TFDs are implemented according to the S-method, as 
described in Article 6.2, the generation of cross-terms can be avoided [8]. 

5.4.7 Summary and Conclusions 

A PWVD of degree p is derived from a p th -order CFD approximation to the deriva- 
tive of the instantaneous phase. It reduces to the ordinary WVD if p = 2. It is real 
and symmetrical about the IF law of an FM signal whose instantaneous phase is a 
polynomial of degree not exceeding p (i.e. whose IF is a polynomial of degree not 
exceeding p— 1), even if that signal is also amplitude-modulated. 

This topic is developed further in Article 5.5 (next). More properties of PWVDs 
are given in [1], Some implementation issues are discussed in [4] and [7]. The effect 
of noise is considered in Article 10.4. 
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